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The variable range hopping results for noninteracting electrons of Mott and Shklovskii are gen- 
eralized to ID disordered charge density waves and Luttinger liquids using an instanton approach. 
Following a recent paper by Nattermann, Giamarchi and Le Doussal [Phys. Rev. Lett. 91, 56603 
(2003)] we calculate the quantum creep of charges at zero temperature and the linear conductivity 
at finite temperatures for these systems. The hopping conductivity for the short range interacting 
electrons acquires the same form as for noninteracting particles if the one-particle density of states 
is replaced by the compressibility. In the present paper we extend the calculation to dissipative 
systems and give a discussion of the physics after the particles materialize behind the tunneling 
barrier. It turns out that dissipation is crucial for tunneling to happen. Contrary to pure systems 
the new metastable state does not propagate through the system but is restricted to a region of the 
size of the tunneling region. This corresponds to the hopping of an integer number of charges over a 
finite distance. A global current results only if tunneling events fill the whole sample. We argue that 
rare events of extra low tunneling probability are not relevant for realistic systems of finite length. 
Finally we show that an additional Coulomb interaction only leads to small logarithmic corrections. 

PACS numbers: 72.10.-d, 71.10.Pm, 71.45.Lr 
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I. INTRODUCTION 

In one dimensional systems, the interplay of interac- 
tion and disorder gives rise to a variety of interesting 
physical phenomena^iL^di^^iMj^. In and below two di- 
mensions, disorder leads to locaUzation of all electronic 
states^iiSiii^iiiii^ii^. Intcractlons lead to a breakdown of 
the Fermi liquid concept in one dimension and to the 
emergence of a Luttinger liquid insteadl'"*'^-, which is 
characterized by collective excitations. In Id electron 
systems, these strong perturbations compete with each 
other. In the presence of a repulsive interaction, dis- 
order is a relevant perturbation and drives the system 
into a localized phase^;^. As a result the linear con- 
ductivity vanishes at zero temperature. Thermal fluc- 
tuations destroy the localized phas o^^'^^ . At sufBciently 
high temperatures, i.e. if the thermal de Broglie wave- 
length of the collective excitations becomes shorter than 
the zero temperature localization length, the linear con- 
ductivity exhibits a power law dependence on temper- 
ature with interaction dependent exponents^. At very 
low temperatures noninteracting electrons show hopping 
conductivityi^ii^ii^. In the present paper we study both, 
the nonlinear conductivity due to strong electric fields 
in the zero temperature localized phase as well as the 
finite temperature hopping conductivity for interacting 
systems. 

Interaction effects in Id electron systems can be de- 
scribed by the method of bosonization^*', which directly 
describes the collective degrees of freedom associated 
with the oscillations of a string. Due to the quantum 
nature of the system, the string can be viewed as a 
two dimensional elastic manifold in a three dimensional 
embedding space. In the localized phase, this mani- 



fold is pinned by the impurity potential. For a classi- 
cal pinned manifold, thermally activated hopping over 
barriers leads to thermal creep, i.e. a nonlinear current- 
voltage characteristio^iiS^. Due to the quantum nature 
of the bosonized electron system, thermal hopping is 
replaced by quantum tunneling giving rise to quantum 
creep2^. 

From a quantum mechanical point of view, tunneling 
driven by a potential gradient corresponds to the decay 
of a metastable state. Using a field theoretical language, 
this decay is described by instantons^i^i^i^l, i.e. so- 
lutions of the Euclidean saddle point equations with a 
finite action. In the presence of a periodic potential, 
instantons relate ground states which differ by a mul- 
tiple of 27r. In this way, the nucleation of a new ground 
state takes place. The instanton picture was used by 
Maki'^^ to study the nonlinear conductivity of the quan- 
tum sine Gordon model. Recently, this framework was 
used to describe quantum creep in a sine Gordon model 
with random phase shifts^-, which models the physics of 
disordered and short range interacting Id electrons in the 
localized phase. 

In the quantum sine Gordon model, instanton forma- 
tion is followed by materialization of a kink particle. The 
materialization is characterized by a rapid expansion of 
the instanton in both space and time. During this expan- 
sion process of the new ground state, the kink particle 
moves according to its equation of motion^^. Within the 
Euclidean action approach, the instanton expansion is de- 
scribed by an unstable Gaussian integral over the saddle 
point, which gives rise to an imaginary part of the parti- 
tion function describing the instability of the metastable 
state^^. In the disordered case, instanton formation de- 
scribes tunneling from an occupied site close to the Fermi 
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surface to a free site close to it. However, spatial motion 
of the kink particle after its materialization is impeded 
by new barriers. 

The central point of our paper is a discussion of the 
materialization process in the disordered case: an exten- 
sion of the instanton is only possible in time direction due 
to the presence of barriers in space direction. The inclu- 
sion of dissipation22i^i^i^iM is essential for this process 
to occur. A dissipative term in the action is essential for 
finding a saddle point associated with an unstable Gaus- 
sian integral on the one hand. On the other hand, dis- 
sipation is necessary to correct the energy balance. It is 
this aspect which was missing in the earlier publication^^. 
The external electric field provides for the energy neces- 
sary to tunnel from an occupied site below the Fermi 
level to an unoccupied site slightly above it. Generally, 
the gain in field energy will be larger than the energy 
difference of the two levels. As dissipation corresponds 
to the possibility of inelastic processes, the excess energy 
can be absorbed by the dissipative bath. In the absence 
of an additional bath, the conductivity was suggested to 
decay in a double exponential manner with decreasing 
temperature^. 

The discussion of the instanton mechanism is combined 
with an RG analysis of the disordered systemi-^-. As 
a consequence, only effective parameters obtained from 
RG enter the results. The combination with an RG is 
necessary to scale the system into a regime where the 
quasiclassical instanton analysis is applicable. 

We make a detailed comparison to MotiM^ variable 
range hopping of noninteracting particles and generalize 
an argument due to Shklovskii'^^ for tunneling of localized 
electrons in strong fields to arbitrary dimensions. The 
instanton mechanism discussed above describes the same 
type of physics. We also discuss the influence of long 
range interactions^. 

We compare the physics of quantum creep to a 
field induced delocalization transition discovered by 
Prigodin'^-i^Si^i^i^. This transition is due to the en- 
ergy dependence of the backscattering probability from 
a single impurity. As the energy depends linearly on po- 
sition in the presence of an external electric field, the 
localization length acquires a spatial dependence. If the 
localization length grows sufficiently fast, the wave func- 
tion is no longer normalizable and the state is delocal- 
ized. We show that this mechanism becomes effective 
at field strengths which are parametrically larger than 
the crossover field for which quantum creep becomes im- 
portant. In addition, we discuss a generalization of this 
mechanism to interacting systems. We find that for re- 
pulsive short range interactions the delocalizing effect of 
an external field is less pronounced (in our approximation 
it even completely disappears) and hence delocalization 
is not a relevant competition for quantum creep. For at- 
tractive short range interactions, the delocalizing effect 
of an external field is enhanced. 

The organization of the paper is as follows: in section 
II we set up the model and discuss its renormalization. 



The tunneling due to instantons and the subsequent time 
evolution of the quantum sine Gordon system is discussed 
as a reference point in section III. Based on this analysis, 
tunneling in the disordered case and the essential influ- 
ence of dissipation is presented in section IV. In section 

V we compare with previous results for noninteracting 
systems obtained by Shklovskii and Mott, in section VI 
we study the influence of a long range Coulomb inter- 
action, and in section VII we discuss the relevance of a 
localization-delocalization transition in external fields to 
the problem of quantum creep in Id systems. 

II. MODEL AND ZERO-FIELD 
RENORMALIZATION 

A. The model 

We consider in this paper disordered charge- or spin 
density waves or Luttinger liquids with a density^'^'^S 

P= (^Po + ^d^^p^ {l + Acos{pip + Qx)) . (2.1) 

Here p — 1,2 for charge density waves (CDW) and Lut- 
tinger liquids (LL), respectively, and ppo = {1/tt)Q. 
The form (|2.ip preserves the conservation of the total 
charge under an arbitrary deformation <p(a;), provided 
|V(/3| ^ Q. The wave vector of the density modulation 
Q = 2kF for LLs and CDWs (although the value of Q can 
be different for CDWs depending on the precise shape of 
the 3D Fermi surface and the optimal electron phonon 
couphng), and Q = Akp for spin density waves (SDWs). 
A = 2 for LLs with a strictly linear dispersion relation^S 
whereas for CDWs y/po^ denotes the amplitude of the 
order parameter of the condensate. 

The Hamiltonian of the corresponding disordered sys- 
tem is then given by 

+2np{x){-Exeo + UR{x))Y (2-2) 

where II(a;) denotes the momentum operator conjugate 
to —ip{x): 

[il{x),ip{x')] :^iS{x^x'), (2.3) 

V denotes the velocity of the phason excitations, and K is 
the measure of short range interactions. The quantities 
can be expressed as = vK — ^^^^ and hence 

_ T h pan ^-^^ „2 _ _po_ .^^]^gj.g ^ jg tj^g compressibility 
and m the effective mass of the charge carriers (for more 
details see Appendix |^ . 

For noninteracting spinless fermions K ^ 1 and v ^ vp 
where vp denotes the Fermi velocity. In CDWs K is 
small, of the order 10~^ whereas in SDWs K may reach 
one^. E denotes the external electric field. 
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The random potential Ub.{x) is considered to be build 
of isolated impurities at random positions Xi, i — 



1, 



1=1 



(2.4) 



In the following, we will neglect the forward scattering 
term ~ dj;ipUji{x) since it does not affect the time de- 
pendent properties. Since disorder is assumed to be weak 
and hence effective only on length scales much larger than 
the impurity spacing, we may rewrite the backward scat- 
tering term in the continuum manner: 

—poAuQ / dx S{x — Xi) cos{pip{x) + Qx) 

i— 1 

-poAuonVpfl^^/^ / dxcos {pip + a{x)) . (2.5) 



Here nimp = Ni^p/L is the impurity concentration, and 
a{x) is a random phase equally distributed in the interval 
< a < 2tt and e*("(^)-"(2')) = aS{x - y). For weak 
pinning, the replacement (j2.5p is valid on length scales 
much larger than the mean impurity spacing and can be 
justified by the fact that both parts of (|2.5p lead to the 
same pair correlations , e.g. the same replica Hamiltonian 



-{poAuo)'^ni 



dxcoa{p{ipa{x) - (pp{x)) . (2.6) 



and hence to the same physics. 

The effective Hamiltonian can therefore be rewritten 



— 2eQLpE — 27rpoAnVpa ^^^uo cos (j)Lp -f a{x)) 

(2.7) 

Below we will use an imaginary time path integral for- 
mulation. It is convenient to use dimensionless space and 
time coordinates by changing x/a — s- x and vr/a y 
where a is a small scale cut-off. The Euclidean action of 
our system is then given by 



S 

h 



K/2T 

dx j dy[{dy^f + {d^^f 

-K/2T 



2t:K 

— 2 f ip — 2u cos {ptp -\- a{x))^ . (2. 
The dimensionless parameters of the theory are 



Ka , ^ KeoEa"^ 

hpv nv 

^ 27ray^KpoAnliluo 

hv 



(2.9) 



(3 is the inverse temperature. All p-dependent physical 
quantities Q{p] K, T, u, /) can be calculated from the case 
p — 1 by the relation 

Q{p;K,T,uJ) = Q{l;p^K,p^T,p\pf). (2.10) 

B. Renormalization group analysis 

Our strategy to consider the transport properties of 
the present one-dimensional system includes two steps. 

(i) We first integrate out phase fluctuations on length 
scales from the initial microscopic cut-off a = 1 to ^ = e' 
1 < < ^ keeping / = 0; ^ is determined from the 
condition that the renormalized and rescaled coefficient 
u{l) of the nonlinear term in (|2.8p becomes of order of 
one or larger and, typically, K becomes small. 

(ii) In the second step we treat the problem at nonzero 
/ in the quasi-classical limit. Corrections to this quasi- 
classical limit arise if the renormalized K is still of order 
one. 

Since the first step is well documented in the literature 
we will here quote only the results. At zero temperature 
and / = the flow equations are given by^iil 



dK 
If 

Boix) 



-^p^u^KBoip^K) 



drr'^e-^ 1 



2\ -^/4 



du^ ( \ o 

dl \ 2 



(2.11) 



(2.12) 



For completeness we also quote the flow equation in 
the case when the phase a{x) is fixed at a = 0: 



dK _ 1 



p^u'B^ip'K), 



(2.13) 



B2{x) = I dr dx {x^ + T^)e ^(1 + y) cosx. 



-x/4 



(2.14) 





-dl - 



Note that in this pure case the definition of u in 
contains an additional factor ^nimpO- 

In both cases (pHl \TT^ and [SUD, there is a 

phase transition a,t K = Kc{u) with Kc{0) = 6/p^ and 
Kc{0) = S/p'^, respectively. For K > Kc the potential 
becomes irrelevant and the system is in a superconduct- 
ing phase. For K < Kc the asymptotic i?G-flow is to 
small values of K and large values of u. For the rest of 
the paper we will restrict ourselves to the case K < Kc- 

The renormalization was done so far for zero external 
field. Now we switch on the field which is a relevant per- 
turbation and destabilizes the system. As in the previ- 
ous section we consider on scales L < ^ only fiuctuations 
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inside one potential valley. The field is then rescaled ac- 
cording to 



di 



2/(0 



(2.15) 



According to our strategy we stop the flow at I = I*, 
e'* = t p'^u{l*) = M with > 1 (see Fig. [l]). Clearly 
our one loop flow equations break down as soon as u{l) 
1, but qualitatively we expect that the flow continues to 
go into the direction of large u and simultaneously small 
K such that typically K{1*) <^ 1. The restricted validity 
of the one-loop flow equations will then merely result in 
the uncertainty of the definition of the value K{l*{Ai)). 




K(l*)=Keff K(0) 



FIG. 1: Flow diagram m K — u plane. RG is stopped at scale 
I = /* such that p^u{r) = > 1. 

The quantities u(l) and K{1) are the renormalized and 
rescaled parameters. The corresponding effective pa- 
rameters observed on scale L = e' are unrescaled. In 
pure systems Kcs{L) = K{1), Ucs{L) = e~^'w(0, and 
/cff(i) = / = e-'^^f{l). In this way we get Ucs{£.) ~ 

For impure systems the corresponding relations are 
ifeff(i) = K{1), ulsiL) = e-^'u^il), Us{L) = / - 
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f{l) and we get u^fj(0 ~ M/{p Sr). Not too close 



to Kci (, can be written as 



FL 



(2.16) 



where LpL = a(p^it)^^/'^ denotes the Fukuyama-Lee 
length. ^ has the physical meaning of a correlation 
length. Note, that in the disordered case there is no 
renormalization of 



V_ _ ^^cff _ 1 

because of a statistical tilt symmetry^. 



(2.17) 



III. TUNNELING - IN THE PURE CASE 

At low but finite temperature the lifetime F^^ of a 
metastable state is determined by quantum-mechanical 
tunneling and given by the relation 



-2ImF = 2f3~^lmlnZ , 



(3.1) 



where F is the free energy. In the path integral formula- 
tion Z is given by 



Z = Zq + iZi 



-i / dyc{{v}J) 
Vlp{x, t) e -^^eff/2T 



(3.2) 

Here /) denotes the Lagrangian of the metastable 

system from which we integrated out already fluctuations 
on length scales smaller than ^. This resulted in the re- 
placement u — > Ucff(^*) and K — s- Kcfi{l*). f represents 
the driving force. The contribution Zq of the partition 
function Z contains the contribution from field config- 
urations restricted to stay in the vicinity of the poten- 
tial minimum. Configurations which leave the metastable 
state are unstable and contribute to the imaginary part 
of the free energy, which is assumed to be small such that 



F « 2im' 



-1^ 

Zn 



(3.3) 



What we will do in the following is the calculation of the 
tunneling probability in the quasi- classical approxima- 
tion: assuming that K^s — K{1*) ^ 1 we assume that 
the only deviation from a completely classical behavior 
is the generation of instantons which trigger the decay 
of the metastable state. Before we come to the disor- 
dered case we briefly exemplify the physics on the pure 
sine-Gordon model which facilitates the further discus- 
sion. The imaginary time Lagrangian of the sine-Gordon 
model is given by 



Lai a 



dx- 



1 



dip 
dy 



dip 
dx 



-2wcff cospip — 2fp 



(3.4) 



The decay of a metastable state (say po = 0) of the sine- 
Gordon model has been first considered by MakiS^ follow- 
ing earlier considerations by t'Hooft^^ and Coleman^-. 
In the present framework the tunneling process is repre- 
sented by the formation of a two-dimensional instanton 
which obeys the Euclidean field equation 

-^p + - pucssmpp + f = (3.5) 

and has a finite Euclidean action. The instanton solution 
Ps{x,y) we are looking for forms a (spherical) droplet in 
the new metastable minimum pi = 27r. Using spherical 
coordinates, = x'^ -\- y'^, we get from p.5p 



(3.6) 



"5^"^+ --^f ~ PUcS smpp -\- f = 0. 
ar^ r or 



For small / and hence large droplet radii R we may drop 
the second term. In particular for / = the solution is 



p{r) 



■ arctanexp [(i? — r) irueffp'^] (3.7) 
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which connects the original state (p{r) =0 for r — R^w 
with the new (hkewise metastable) state ip{r) ~ 2tt/p for 
R — r ^ w where w = \/ ^Jp^u^s ~ ^/a/AT is the width 
of the droplet wall. The neglect of the second term in 
p.6p is justified for R:^ w. In the following, we will use 
this narrow wall approximation to describe the instanton 
only by its radius ignoring for the moment the shape 
fluctuations. The action of the instanton includes then a 
surface and a bulk contribution 



S{R) 



where a 



2tt 
pKcS 

Pi 



{Ra - R^f) 



(3.8) 



is the surface tension of the 



instanton (the Euclidean action associated with the unit 
length of the instanton boundary is equal to ha / (jjK)) . 
The action (13.81) has a maximum at 



Rc = 'y/{2f). 



(3.9) 



where S{Rc)/h = ■S'inst/^ = t^ct^ /C^pKcsf) denotes the 
instanton action. 

The ratio iZi/Zo appearing in the p.3p can now be 
calculated as a functional integral over circular droplets^ 



^1 T 

— cx Im 



'S'inst 



Here we have taken into account that the change from the 
integration over ip to R involves a Jacobian ~ S-mst/h. 
The integral, taken along the real axis, is divergent since 
the saddle point (in the functional space) is a maximum 
of the action as a function of the instanton size R = 
Rc + r. 

This leads - after the proper analytic continuation - to 
the imaginary part in iZi in the partition function (for 
a more detailed discussion see Schulman^^). As a result 
we obtain 



iZ 



.R 



Zo ? 



(3.10) 



Further contributions to the pre-exponential term result 
from the inclusion of fluctuations of the shape of the in- 
stanton. Within the narrow wall approximation these 
fluctuations could be taken into account by rewriting the 
Euclidean action in the form^ 



5; 



K^a/2T 
-K,tt/2T 



'^l + {dyXkf - 2kfXk{y) 



(3.11) 

Xk{y) (where k — ±.) denote the wall position of the 
left and right segment of the wall of the critical droplet, 
respectively. These fluctuations were considered by 



Maki^Sj they lead to a downward renormalization of the 
surface tension in (|3.1ip . 

So far we considered the center of the instanton to be 
fixed at a; = Xq and y — 0. However the position of 
the center can be moved around which corresponds to 
the existence of two modes with zero eigenvalue. Cor- 
respondingly, in the calculation of the partition function 
we have to integrate over all possible instanton positions 
which delivers in the low temperature limit 



^ R - 



an additional factor^^ 



inst 



2nh 



(3.12) 



(3.13) 



where Lq is the sample length. The first factor counts the 
number of different instanton positions, the second one 
is the Jacobian resulting from the transition from the 
original to the displacement degrees of freedom. Thus 



-Tza''/(2pK,!lf) 



(3.14) 



represents the time scale on which the phase field 
ip tunnels through the energy barrier between tpo and ipi . 

Once the critical droplet of the new ground state ipi 
is formed, the field materializes and evolves subsequently 
according to the classical equation of motional. 

The latter follows from Eq. (|3.5p by the resubstitution 
y — ivt. Since ipsi'r) depends only on = — v^t^, the 
solution of Eq. p.6p also describes the evolution of the 
phase field after the tunneling event. Within the narrow 
wall approximation, the Minkowski action is 



Sm 




fe=± 



Rc J 



(3-15) 

which leads to the equation of motion of two relativistic 
particles 



± 



1 d 



x±it) 



dt 



1-Xl/v^ 



R-\ 



(3.16) 



This has to be solved with the initial condition X±{t — 
0) = Xq ± Rc where t = corresponds to the moment of 
the materialization, hence 



X±{t)^±jRl+vl^f' + XQ 



(3.17) 



Clearly, with ivt — > y, this is also the saddle-point equa- 
tion of action (|3.1ip . corresponding again to a spherical 
droplet of radius Rc ■ 

The energy £ is a conserved quantity during the motion 
of the kinks as 



5 = 



= 0. 



(3.18) 
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Since the potential energy of the new metastable state 
is always lower than that of the initial configuration the 
kinks will accelerate, approaching eventually the phason 
velocity Weff. 

In a long sample the nucleus of the new metastable 

(n) 

state may form independently at several places Xq fol- 
lowed by a rapid expansion of the kinks X^\t) which 
finally merge. Below we will show that the correspond- 
ing picture in a disordered sample is rather different. 

For the sake of completeness and since we will use it 
later, we consider the influence of an additional damping 
term in the Euclidean action^^iSiiH 

-K^!!/2T y T ' L ifrff J 

(3.19) 

We have introduced a phenomenological dimensionless 
constant rj describing dissipation. This form of the Eu- 
clidean action corresponds to a linear (Ohmic) dissipa- 
tion term in the classical limit in real time. 

In the saddle point equation (|3.5p . (|3.19p results in an 
additional term 



2tt 



-K,h/2T 



sm 



2 \ T-K(y-y') l 



(3.20) 



In appendix B, we justify the choice of this dissipative 
term and explain its physical origin. 



Here, a[Xk,dyXk) denotes the surface tension of the in- 
stanton which depends in general both on x and on the 
slope dyXk of the surface element. We could now pro- 
ceed in looking for solutions of the saddle point equation 
belonging to (|4.2p . Since this equation depends on the 
specific disorder configuration a general solution seems 
to be impossible. As will be shown below, the surface 
tension (j[x, O) = ax{x) changes on the length scale ^ by 
an amount of the order cr whereas a(^x,oo) = ay ^ a 
is constant. Changing the slope of the surface element 
we expect an monotonous change of a(^x,dyX) between 
ax{x) and ay. Because of the rapid change of ax{x) the 
number of saddle point solutions will be huge. But only 
one of them will be relevant for the present problem in 
the sense that it rules the decay of the metastable state. 

Following previous workli^i?., we calculate next the sur- 
face tensions ax and ay in the limit Kcft <C 1. This is in 
agreement with our general strategy to take into account 
quantum effects (i.e. the fact that KcS 7^ 0) only when 
instanton formation is considered. Additional shape fluc- 
tuations of instantons - which are also of quantum origin 
- could be included later on by considering the renormal- 
ization of this surface tension. In general they will reduce 
the surface tension by a finite amount as long as we are 
below Kc- As it was shown irtiiii^ the limit K^s ^ 1 
allows an exact solution for the (classical) ground state. 
To reach this goal we rewrite the effective action as a dis- 
crete model on a lattice with grid size ^. In the classical 
ground state, ip{x, y) does not depend on y any more, and 
the renormalized form of Eqs. (|2.71 12. 8p can be written 
as 



IV. TUNNELING - THE DISORDERED CASE 

A. Surface tension of instantons 

In this section we will consider the tunneling process in 
the disordered case. In principle we follow the calculation 
of Section IIIIl but important differences apply. Starting 
point is again the effective action of the form (|2.8p on 



the scale ^ 



ae 



where u is replaced by u(l*) and K is 



replaced by if (r) as follows from (pUj) and ^J^. The 
saddle point equation of the instanton now reads 

^ + ^) ^ - "(^*) ^i'^ (P'fi + "(^)) + /(^*) = • 

^ , (4.1) 

Let us assume that we have found its narrow wall solu- 
tion which is parametrized by the instanton shape X± (y) . 
Plugging this solution into the effective action the latter 
will take the general form 



n 



dy la{Xk, dyXk) ^l + {dyXkf 



K^tt/2T 



-2kfXk{y) 

(4.2 



S_ 
h 



Hctf 



2ttT 



dx ^ [dx'^f - cos {pLp + a{x)) 



(4.3) 

Next, we replace the integration over a; by a summation 
over discrete lattice sites i = x/£,. This gives 



Hoif 1 



T 2ttT^ 



{{(pi-(pi+i)'^-2u{r)cos{p(pi+2TTai)} , 

i=l 

(4.4) 

with < cti < 1, and 2TTai denoting a random phase. 
If u{r*) = Ai/p^ 00, the exact ground state can be 
written afM^ 

p(fi{m) — —2'Kai + 27r^^[Q;j — Ofj^ijc + 2'Km (4.5) 

where \o\g denotes the closest integer to a, and m is 
an integer. A uniform electron density corresponds to 
a constant (pi(rri). Disorder leads to an inhomogeneous 
electron distribution with the maximum of | (^i+i — (pi \ 
equal to 7r/p corresponding to a maximal excess charge 
of ±eo/2 on the scale ^. 

Next, we consider a possible bifurcation of the ground 
state. Indeed, since the (fi are constructed such that local 
energy is always minimized in the ground state, we have 
merely to minimize the elastic energy, which depends on 
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the square of 

— {<fi{m + n)-(fi^i{m)) = - {5a^ - [dai]G) + n, (4.6) 

where n is integer. Since Sai = OLi — cti-i is equally 
distributed in the interval —2 < 8ai < 2, we have 
-1/2 < Sa, - [6a,]G < 1/2. For \Sa^ ~ [Sa^U < 1/2 
the ground state is uniquely determined by n = 0. How- 
ever, for Sat — [Sailc = k/2, k = ±1, phase config- 
urations with n = and n — k have the same en- 
ergy and hence the ground state bifurcates. The pairs 
of sites which show this property can take the values 
ifi — (fii-i = ^kn/p. Bifurcation corresponds to a lo- 
cal change of the charge by ±2eo/p. Since adding (or 
removing) a charge (or a pair of charges in the case of 
CDWs) does not change the energy one has to conclude 
that these bifurcation points correspond to states at the 
Fermi surface. The states which fulfill exactly the bifur- 
cation condition have measure zero, but there is a finite 

< e. For those pairs the cost 

for a deviation from the ground state is also of the order 
e. 

To calculate the surface tension we first consider a wall 
parallel to the y-axis. Such a wall - at which m is changed 
to m + n, n = ±1 - has an excess energy 

-(^,(™)-^,_i(™))'}= (4.7) 



probability that 



On the rhs we introduced the surface tension of the ori- 
ented (n = ±1) wall Uxi(ji) = avi(n), a « 2i:/(jpS) 
and2^ 



9i = (P, ( X i ) 




2jt/p 




-EeQX 



FIG. 2: Upper figure: two equivalent ground states witfi two 
bifurcation sites which facilitate the formation of an instanton 
in nonzero external field. Middle: Top view of an almost 
rectangular instanton, the width of the wall corresponds to 5, 
i.e. the distance between the grid points of the upper figure. 
Lower figure: Hopping process corresponding to the instanton 
formation, see also Fig. 3 in^. 



v,{n) ^l~2n{5a,-[5a,]G) , + ^ 2 (4.8) 

and hence < i'i{n) < 2. Thus we find that if we pro- 
ceed in cc-direction the surface tension ax,i{n) of a wall 
with a fixed orientation is fluctuating from site to site 
and equally distributed in the interval < ax.iin) < 2a. 
If the surface tension for a -I— wall is close to 2a the cor- 
responding surface tension for the minus wall is close to 
zero. If we consider the surface tension of an oriented 
wall in a typical region (i.e. we exclude for the moment 
regions where rare events take place) of linear extension 
Lx 3> then the average m-th lowest value of i'±{x) in 
this region is given by 2171 

An analogous calculation can be done for the surface 
tension ay if we introduce also a lattice of grid size ^ 
in the j/-direction. Since the disorder is frozen in time, 
Sai = and hence ay{n) — a. Note that the value of 
a found here agrees up to a numerical constant with the 
surface tension obtained for the pure case in the previous 
section. 



B. Instanton action in the disordered case 

To describe the decay of the metastable state - which 
is in the present case one of the classical ground states 
(pi{m) equation ()4.5p - the functional integration in Zi 
has to include the integration over droplet configurations 
of the new phase 'f>i{m -t- 1) starting with small droplets 
which subsequently enlarge in both x and y-direction un- 
til the new metastable state is reached at least in a part 
of the \d system. Since the surface tension ax can be- 
come arbitrarily low it is obvious that the dominating 
droplet configurations are those bounded by walls essen- 
tially parallel to the y-axis. Once the saddle point (the 
instanton) is reached, the droplet will expand only in 
y-direction since further expansion in x-direction is pro- 
hibited by the high cost of leaving the wall position with 
extra low surface tension. The maximal extension of the 
droplet in a;-direction will be denoted by Lx^c- 

We will start our quantitative consideration by restrict- 
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ing ourselves to rectangular droplets of linear extension 
and Ly, respectively. The instanton action S{Lx^ Ly) 
is then given by 



) pKcfi 



S{L,,Ly)^Lx{2 + r]ln^}- 



h a 

(4.9) 

where we have included also the contribution of the 
dissipation^. We will assume here that the dissipation is 
weak, T] It should be noted that the dissipation fa- 

vors in general a rectangular shape of the instanton even 
in the absence of disorder—. v±{x) denotes the gener- 
alization of = ±1) to a continuum description and 
is a random quantity equally distributed in the interval 
< < 2 with a correlation length of the order ^. is 
still given by ()3.9p . The functional integral corresponding 
to p.lOp takes now the form 



(4.10) 



Note, that the second equation describes the energy con- 
servation during tunneling. These equations have the 
solution 



l + §ln(^) 

l + iMf) 
ri ( eR, 



- i?c ( 1 + ^ln^ ^ 



Thus we get for the saddle point action 
Sr 2a ' 



(4.13) 



(4-14) 

This solution becomes exact in the quasiclassical limit 
Kcff ^ 1. For < 1, there are additional fluctuations 
which lead to downward renormalization of the surface 
tension. 



C. Stability of the saddle point solution 



It is illustrative to consider the evaluation of (|4.10p for 
the pure case, where v±{x) = 1, and ry = 0, ignoring for 
the moment the fact that the consideration of circular 
droplets is here more natural. The main contribution 
to (|4.10p comes now from quadratic droplets since for a 
given area the square has the least circumference. At the 



saddle point L^., 



2Rc and the instanton action 



is Sc — 4i?c- Clearly, circular droplets have a lower saddle 
point action Sc = tt Re- 
in the impure case the situation is different because of 
the variation of the surface tension i'±{x). The saddle 
point will be determined by low values of the surface 
tension v+{x) + v-{x + Lx). Since the surface tension 
v{x) fluctuates on the scale ^, the droplet will not - as 
in the pure case - expand both in the x- and in the y- 
direction. For the further discussion it is useful to use 
the parametrization — A ■ where < < 1. In 
an x-interval of extension A the typical surface tension 
fulfills the inequality 

2^<iy+{x)+iy.ix + Lx) = 2^fi{x,lx)<2. (4.11) 

For a typical region ^{x,lx) ~ fJ-{lx) will not depend on 
the choice of x. As explained above ii{lx) is strongly 
fluctuating with a correlation length ^/A. Since we have 
no detailed information about the strongly fluctuating 
surface tension fi{lx) we now minimize S, Eq. (j4.9p . with 
respect to the scale factor A and Ly, and choosing such 
a value of lx that fi{lx) is of the order one. This gives the 
saddle point equations: 



Lx{2 + 77ln ^) - {2^fi{lx) + ^)Ly = 0, 



L, 



11- 



2^m(^x) 



Lx 
Lin 



(4.12) 



0. 



Considering quadratic fluctuations around the saddle 
point it is easy to show that one eigenvalue is negative as 
it has to be the case for an instanton. The eigenvector 
corresponding to this eigenvalue has a small component 
in x-direction which could suggest a further growth of 
the droplet both in y- and in x-direction. For the fur- 
ther discussion it is however important to realize that 
the decay from the true saddle point has to follow one of 
the many surface tension minima inscribed in the func- 
tion fi{lx) and which determines the saddle point value 
of lx- The latter would follow from a more microscopic 
approach using a variation with respect to Lx = Alx- 
Such a treatment would require the detailed knowledge 
of the energy landscape fJ.{lx) which is unknown. Since 
fi{lx) is of the order of one for the low surface tension 
valleys we conclude the for the true saddle point Ix^c = Ic 
and fi{lc) = IJ-c remain unspecified, but are both of order 
unity. Such a decay within one energy valley can only 
happen if (i) the slope of S in y-direction is negative and 
(ii) if 



dLl 



I l^-u — -^'i; 



(4.15) 



The first condition is fulfilled, as one can see from the 
saddle point equation for Ly. The second condition de- 
pends crucially on the existence of dissipation. Without 
dissipation, e.g. the emission of phonons, the decay to 
the new metastable state is impossible because it would 
violate energy conservation. In experimental systems, 
dissipation is always present, for example due to electron 
phonon coupling. The precise way in which the dissi- 
pation strength enters the final result as a prefactor is 
beyond the logarithmic accuracy of our present calcula- 
tion. However, from a comparison with the decay rate of 
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dissipative two level systems^ we conjecture that the de- 
cay rate is proportional to the damping strength rj. The 
answer to the question whether an short range interact- 
ing and disordered ID electron system has a sufficient 
amount of intrinsic dissipation to support a finite creep 
current is clearly beyond the scope of the present work. 

As a final result we get for the decay rate T{E) in the 
limit of 77 <C 1, E <^ Eq and hence for the nonlinear zero 
temperature conductivity <j{E) 



\nT{E)^\na{E)^^2^^ fl + r;ln^^^^^ 



where 



En 



doc 



(4.16) 



(P'^eoCioc) 



Cioc V E pKea^f^^ 

(4.17) 

Lx{E) denotes the extension of the instanton and hence 
the distance over which charges are tunneling in x- 
direction. k is the compressibility which depends on 
the interaction parameter K (compare Appendix A). 
60c = P^£,KeS is the localization length of the tunnel- 
ing charges, in agreement with a result of Fogler— . In 
fact the correlation length ^ was used in^-^a a^d subse- 
quent papers as the localization length, which is correct 
for intermediate K. In the limit K 0, ^ goes over into 
the Fukuyama-Lee length which is a completely classi- 
cal quantity. Tunneling processes are characterized by 
a length scale ^oc which is intrinsically of quantum me- 
chanical origin and therefore has to vanish for ft — > 0. 



D. Percolating instantons 



Once a single instanton was formed it will expand in 
the y-direction leaving a region of linear extension (E) 
behind in which the phase is advanced by 2-K/p. Since 
such a configuration corresponds to a change of the phase 
difference of neighboring sites at the left and right bound- 
ary of this region by -izlii /p one has to conclude that such 
an event corresponds to a transfer of a charge 2eo /p from 
one boundary to the other. Many such events will hap- 
pen independently at different places of the sample, each 
of it leads to a local transfers of charge over a typical dis- 
tance Lx{E). A current will only flow if these tunneling 
processes fill eventually the whole sample, see Fig. 3. 

It is in this respect important to note that bifurcation 
sites (at which the original metastable state (^j(m),Vj 
was changed to ipj>i = (pj{m + n) [n = ±1)) are now 
favorable sites for a change to Lpj — (p,j{m + n),yj. In- 
deed, if the initial reduced surface tension for this pair 
of sites was (compare (|4.7p ') i^i = e, then we get for 
the surface tension between the state (fj (m + n) , Vj and 
iPj>t = (pjim + n),(pj^i = (pj{m) 




FIG. 3: Phase profile fix) (bold line) after several tunneling 
events have happened, starting from the zero field ground 
state characterized by the dashed line. The thin lines show 
different zero field ground states which follow from each other 
by phase shifts of 2-K/p. A current flows only if the instantons 
form a percolating path through the sample. 



T 



-(^,(m + n)-^i_i(TO))^| = (4.18) 
The reduced surface tension of the new wall is therefore 



Vi{n) = -1 + 2n((5a,; - [Sai]G) = 



(4.19) 



i.e. these sites are preferential for the formation of new 
instantons which will fill the gaps between already exist- 
ing regions of increased Lp. These considerations have a 
simple interpretation in terms of the charges which un- 
dergo tunneling: In the initial ground state (pi{m) all 
occupied states have energies below the Fermi level. The 
creation of an instanton corresponds to the transfer of a 
charge (cq for Luttinger liquids and 2eo for CDWs) from 
an occupied to an unoccupied site of distance (E) . The 
site which is unoccupied now is below the Fermi-energy 
and consequently its surface tension is negative. 

As one can see from (j4.16p the current / ^ ~'^{'f) ^ 
T{E) shows for E Eq creep-like behavior. The calcu- 
lation of the complete pre-factor of the decay rate F is 
beyond the scope of the present paper. 



E. Finite temperatures 



So far we assumed that the typical extension Ly,c of 
the instanton in y-dircction is much smaller than K^g/T. 
However, this condition is violated for 



f3<(3i 



P 



eoE^ 



E <E, 



loc 



/3 



P 



eo6oc/? 



(4.20) 



For higher temperatures, P < (3e, the saddle point can- 
not be reached any more. In the limit /3 ^ /J^ we can 
proceed as in2^ and consider the contribution to the av- 
erage current j ^ —^{ip) from the production of droplets 



10 



of the metastable states with n = ztl.This leads in the 
hmit PcqEL^ 1 to a hnear relation between j and 
E. The maximal barrier is now given by Ly « Kos/T, 
Lx ~ Lx{Ep) which gives a linear conductivity 



a{T) - exp { - const. ^ftwp/J/Xeff} • (4.21) 

Here, T — una/ [3 where [3 is the inverse temperature, 
and we introduced the pinning frequency Wp by 



The result (|4.2ip can be also written in the form 



cr(T) - exp { - const. \/^/«ioc)}. 



(4.22) 



(4.23) 



where we used the compressibility k defined in the Ap- 
pendix. Note that n depends on the interaction param- 
eter K. Equation (|4.23p is the natural extension of the 
Mott variable range hopping result to short range in- 
teracting electrons. A long range Coulomb interaction 
leads to a logarithmic correction in Eq. (|4.23p . see^. 
The additional logarithmic correction due to a dissipative 
damping of the tunneling action by electronic degrees of 
freedom^ is only relevant for strong pinning, whereas in 
the situation of weak pinning this correction is absent as 
all electronic degrees of freedom are localized and cannot 
dissipate energy. 



F. Rare events 

So far we considered typical instantons. The question 
arises what is the influence of regions in which the lowest 
surface tension is untypically large? The question was 
addressed for nonzero temperatures^i^i^. We proceed 
in a similar manner. The regions with large surface ten- 
sion will lead to a lower tunneling probability and hence 
act as a weak link. We may first ask the question: what is 
the largest distance Xm between two consecutive sites of 
surface tension Vi < m a sample of length Lo? These 
sites are randomly distributed along the sample of length 
Lq. Their concentration is given by 



"0 = ^'o/C • 



(4.24) 



The probability distribution for the separation x between 
such consecutive sites is 



(4.25) 



The number of sites with Vi < has the average LqUq 
and a standard deviation ^/Lquq. Therefore, the largest 
separation Lx^m between such sites can be found from 



Lo"-o / dxpo{x)^l. 



(4.26) 



where the integral represents the probability that the 
given pair of neighboring sites has the separation exceed- 
ing Lx,m- Using Eqs. ([424| and (|425ll . we find: 



-—\niLo- 



In 



L.r 



(4.27) 



Comparing this equation with (|4.1ip we see that in the 
region considered fi{x,lx) is of the order hi{Lo/ Lx,m) in- 
stead of order one as in the typical regions. Combining 
this result with the saddle point equations (|4.12|) we find 
that we have use the replacement 



fi{lx) pi.{lx) ln(£o/A 



(4.28) 



in (j4.13p and the subsequent formulas. Since it is this 
weak link which will control the total current ,the main 
exponential dependence of the current is given by 



In/. 



\ 




(4.29) 



We may now consider the influence of this weakest link 
on the conductivity at finite temperatures. We found 
before that the nonlinear conductivity crosses over to a 
linear behavior when the extension of the instanton in 
y-direction reaches Kf.ff/T = Ly^c ~ Rc- This relation 
is not affected by the value of the surface tension as can 
be seen from (|4.13l) . We have therefore to use the same 
replacement E — ^ Ej^ as for the typical instantons but 
take into account the extra factor in (I4.29p . which gives 
for the resistance of the the sample of length Lq 



HRrn/Ro) 



loc 



/In io 



/36 



loc 



(4.30) 



Finally we mention that for exponentially large sam- 
ples Lq ~ ^6^"=/^ there are weak links with surface tension 
1^0 of the order one. In this case the tunneling probabil- 
ity is \nT{E) ~ Eq/E. A similar crossover with respect 
to the sample length was found in Refs. |46||47| for linear 
conductivity in Id in the case of nonzero temperature. 



V. SHKLOVSKII AND MOTT VARIABLE 
RANGE HOPPING 

In this Section we give a brief account of the deriva- 
tion of Shklovkii's zero temperature nonlinear hopping 
conductivity which was done for three dimensions'^^. Its 
extension to general dimensions d is straightforward. All 
electrons are assumed to be localized over a distance of 
the order of the localization length ^loc- An electron 
can be transferred from a filled to an empty site sepa- 
rated by a distance Lx without absorption of a phonon 
provided the difference between the difference of the en- 
ergies of the states does not exceed e^ELx- If we de- 
note the density of states (per energy and unit length) 
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by g{£), the number of states states with energy smaller 
than cqEL^ accessible by a jump over the spatial dis- 
tance is given by giSp) ■ bqEL'^^^ > 1. £jr denotes 
the Fermi energy. The minimal distance for the tunnel- 
ing is therefore given L^^c ~ (eo£^3(^-F)) The 
probability for an electron to jump over a distance L^^c 
is P{Lx) ~ g-2i^,c/4ioc ^ / aj^j hence we get for the 
nonlinear current voltage relation 



(5.1) 



This result agrees with our findings (|4.16p in d = 1 
dimension in the noninteracting case K — 1. 

Indeed, the one-particle density of states at the Fermi 
level can be written as 



^0 — 



(5.2) 



where are the energies of the eigenstates of the Hamil- 
tonian describing noninteracting electrons, and Lq is the 
system size. At zero temperature, the particle density p 
can be expressed as 



Lo 



e-k) , 



and hence k = dp/ dp. = g{p) for p, = £p. 

In the case of finite temperatures when eo-E6oc/3 < 1 
the number of reachable sites is determined by thermal 
activationi^iii. Indeed, the number of these sites located 
in a volume and the energy interval AE is given by 
g{£F)AEL'^ > 1. The hopping probability is then to be 
determined from the maximum of 



c/6ocg-/3/(s(£F)i^) 



which gives 



J^x,l3 = I 



\2g{£F) 
and hence for the current 



I ^ E ■ exp 



2(d+ 1) 



/ pel \ 
y2g{£F)^iJ 



i/{d+i) 



(5.4) 



(5.5) 



(5.6) 



Thus there is a crossover from a nonlinear current volt- 
age relation at low temperatures (3 > (3e = (eo£'^ioc)~^ 
to a linear relation at higher temperatures. The result 
(|5.6p agrees for K — 1 again with our result (|4.2ip . 



VI. THE INFLUENCE OF COULOMB 
INTERACTION 

Efros and Shklovskii^^ considered the influence of the 
Coulomb interaction on Mott variable range hopping. 



The Coulomb repulsion leads in dimensions d > 1 to 
a suppression of the density of states close to the Fermi 
energy, g{£) ^ \£ — £F\'^^^^ which changes the exponent 
l/(d-|- 1) in ((57T1) and ([?^ to 1/2 in aU dimensions. This 
can be seen as follows: 

To calculate the density of states in a system with 
Coulomb interaction we first consider deviations from the 
ground state. For any pair of localized states close to the 
Fermi surface the transfer of an electron from am occu- 
pied site Yi of energy £i < £f to an empty site Yj of 
energy £j > £p the net change of energy has to be posi- 
tive 



A£ = £j -£^- 



> 0. 



(6.1) 



Here Cs denotes the dielectric constant. Let us now as- 
sume that we consider all energy levels with | £i — £f \< 
£. Apparently, they have to fulfill the inequality (|6.1D 
which gives 



l> 



(6.2) 



Hence the donor concentration n{£) cannot be larger 



(5.3) than n{£) < (j^) from which wc find for the den- 



sity of states gi£) = « i'^-^ef/el'^. Thus there is 

no change in the density of states in c? = 1 dimensions. 
A slightly more refined calculation gives a logarithmic 
correction to the density of states which give logarithmic 
modifications of in the exponents of (|5.ip and (jS.G p ^^'*^^ . 

We consider now the influence of long range Coulomb 
interaction. In the Fourier transformed action 
k^\(Pk\ 

dimensionless prefactor C — 



is replaced by k'^{l+C\n{kF/k))\(pk\'^ where the 



measures 

the relative strength of the Coulomb interaction. Fluc- 
tuations of </? are now reduced leading to a suppression of 
the transition to the delocalized superconducting phase. 
The effect of the Coulomb interaction in the localized 
phase is weak and we have still a flow of our coupling pa- 
rameter u to large values. Following our considerations 
of Section IV we have next to construct the exact ground 
state in the presence of Coulomb interaction. Because of 
its nonlocal character, Eq. (|4.5p is not longer the true 
ground state. However, it is clear that any configuration 
which deviates from the ground state has to increase the 
energy. We consider now the instanton as such a con- 
figuration. The effective instanton action is now given 

by 



ha ^ i ' 



^{{p{x,l^)-C)j- - ^]Ly. 



(6.3) 



Here p denotes now the surface tension in the presence 
of the Coulomb interaction and there is an additional 
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contribution ^ C ^ C from the Coulomb interaction be- 
tween the surfaces. Since we consider deviations from 
the ground state, fi{x,lx) — C has to be positive. It is 
therefore tempting to assume that this difference scales 
as ^{x,lx) — C ~ fJ-cilx) and the results of Section IV ap- 
ply with fi{lx) replaced by fJ-c{lx), in complete agreement 
with the consideration for the case without long-range 
interactions. Thus, the Coulomb interaction does not 
change the power law dependence of the the results (|4.16p 
and (|4.23p on E and /3. There is a logarithmic correc- 
tions to nc{lx) due to long range Coulomb interactions'*^ 
which is beyond the accuracy of the present calculation. 



VII. DELOCALIZATION IN STRONG FIELDS 

For our discussion of quantum creep in Id electronic 
systems the presence of dissipation was essential. A 
possible nonlinear I-V characteristic in noninteracting 
electron systems without dissipation was discussed by 
Prigodin^S and others^^iiSiiilii^. This nonlinearity is due 
to a field dependent delocalization of electronic wave 
functions. To understand whether this mechanism pos- 
sibly interferes with the quantum creep discussed in sec- 
tion IV., we briefly review its derivation and discuss a 
generalization to short range interacting electron systems 
afterwards. 

The disordered potential (2.4) consists of a series of 
delta function scatterers of strength uq placed at irregu- 
lar positions Xi with a mean density riimp- A free electron 
incident with energy £ on one of these scatterers is re- 
flected back with a probability 



1 



(7.1) 



Neglecting localization effects for the moment, the en- 
ergy dependent mean free path is given by £{£) — 
i?~^(£)nr^p. In the presence of an external electric field, 
the electron kinetic energy becomes position dependent 
according to £{x) = £o + eEx, leading to a position de- 
pendent mean free path 



iix) 



2ft2 

1 H k{£o + CqEx) 



1 



^imp 



(7.2) 



for noninteracting electrons. In one-dimensional systems, 
the localization length ^loc is approximately equal to the 
mean free pathi^. The equality between mean free path 
and localization length is valid even in the presence of 
an electric field, as the additional phase shift that the 
electron acquires as it moves along the field is canceled 
out on the return trip counter to the field"^^. Thus, the 
energy dependence of the backscattering probability of 
an individual scatterer gives rise to a position dependent 
localization length ^ioc(a;) ~ £{x). 

Without applied external field, a localized electronic 
state centered around a position xq is characterized by 
an exponentially decaying wave function envelope ^>{x) ~ 



exp(— |a; — xo|/^ioc)- In the presence of an external field, 
the localization length acquires a position dependence. 
The electron kinetic energy changes on a length £o/eoE. 
If the localization length is much smaller than this length 
scale, ^loc <C £o/eoE, the localization length varies slowly 
in space, and the wave function envelope is given by^ 



*(a;) 



— dx-p — ^^T^ 

e -^^0 £ioc(^) , X > xo 



(7.3) 



Using the explicit position dependence of the mean free 
path given in (|7.2p . one finds that the envelope decays 
asymptotically as a power 



1 + 



{x — Xf))eE 

£'o 



'""-imp "0 



(7.4) 



The state is only localized if it is normalizable, i.e. 
integral over the envelope squared 



the 



dx \I'^(x) < CO 



(7.5) 



must be finite. The wave function is only normalizable 
if the envelope decays faster than 1 / y/x or if the electric 
field is weaker than the critical electric field 



Ec = 



(7.6) 



This calculation underestimates the exact result^-* by a 
factor of two. We want to compare this critical field to the 
crossover field Eq = obtained from Eq. for 

Sloe 

the creep current. Using the relation ^loc ~ muin'^ valid 
for noninteracting electrons, the crossover field Eq. (j7.6|) 
can be rewritten as Ec = i^o^F^oc ^ Eq. From this re- 
lation one sees that the instanton mechanism is effective 
already at fields which are parametrically smaller than 
the critical field for delocalization. 

If the point scatterers are not of delta-function type 
but characterized by a finite interaction range instead, 
even an arbitrarily weak external electric field will lead to 
delocalization on long length scales as the backscattering 
probability drops significantly for electron wave vectors 
larger than the inverse potential range. 

How can this argument be generalized to short range 
interacting electrons? If the particle density is spatially 
homogeneous, the kink kinetic energy far from the cen- 
ter of localization is much larger than the Fermi energy 
of the interacting system and the kink should essentially 
behave like a free particle. In this case, the kink wave 
function would stay localized for fields below the thresh- 
old value l|7.6p , and the power law tail of the kink wave 
function could possibly enhance the tunneling current de- 
rived in Sec. IV. Alternatively, the charge density may 
be locally in equilibrium with the electro-chemical po- 
tential. In this situation, both the particle density and 
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the Fermi wave vector are position dependent. Assum- 

2 

— is at least 



hik- 



ing that the free electron relation Ep 
approximately valid, the Fermi wave vector acquires a 
position dependence 



2m 

kp (x) = \l {£o + eEx) 



(7.7) 



Expressing the cutoff as A ^ kp, the localization length 
depends on the Fermi wave vector according to 



Cloc OC kp\kpLFi^)TT^T^7Kn 



{7- 



where the classical Fukuyama-Lee length ipL is indepen- 
dent of kp. For a noninteracting system, K/K^ = 2/3 
and one retrieves the dependence ^loc ^ kp. For general 
K, the position dependence of the localization length is 



^ioc(a;) OC 



2m 



eEx) 



(7.9) 



Accordingly, the wave function envelope decays as a 

—const X i<c-K j for K < 1 

and does not decay to zero at all for K > 1. The ar- 
gument leading to this result contains approximations 
which may be correct only qualitatively. For this rea- 
son, we only draw the conclusion that the delocalizing 
effect of an external field is weakened by repulsive short 
range interactions and enhanced by attractive ones. In 
order to determine whether the critical field strength Ec 
indeed has a discontinuity when varying K across one a 
more sophisticated calculation is necessary. 

In conclusion, the delocalization mechanism studied in 
this paragraph does not seem to interfere with the creep 
current described by formula Eq. (4.25) for the case of 
repulsive short range interactions. 



disorder favors the charge transport: the corresponding 
pure system has an excitation energy gap. 

In the case with disorder, the notion of the materi- 
alization is different from that in the pure case. The 
disordered system is intrinsically inhomogeneous, so that 
the transport occurs as a sequence of local materializa- 
tion events. In the fermion picture these correspond to 
separate hops whose positions and lengths depend on the 
disorder configuration. At zero temperature, the typical 
hop length depends on the electric field E and is given 
by Lx = i,ioc\/ Eq/ E. In this regime, the conductivity 
reads: a cx e:cp{~ Eq / E) . Quantum creep takes place 
for E ^ Eq. The results can be expressed in terms of 
the localization length and the compressibility. 

We compared the electric field Eq with the thresh- 
old field Ec found in Refs. [HilliOlilllii: Eq ~ 
Ec/{kp£,\oc) ^ Ec, therefore we expect that our creep 
conductivity result is valid. 

We demonstrated that the inclusion of the dissipative 
term in the action is necessary to find an unstable mode 
in the functional integral. Besides, the dissipative bath 
takes into account the possibility of inelastic processes 
and ensures the energy conservation. The question about 
the energy balance does not appear in the treatment of 
the pure system where the gain in the electric field energy 
is spent on the increase of the kinetic energy of departing 
kinks. 

However, for weak dissipation, our result for the ex- 
ponent in the quantum creep regime only weakly de- 
pends on the dimensionless dissipative coefficient 77. This 
feature can certainly change as one crosses over to the 
nonzero temperature regime with linear conductivity. 

The effect of the Coulomb interaction in ID was shown 
to be irrelevant. It destroys the delocalization transition, 
but it is not expected to change much in the phase al- 
ready localized by the disorder. 



VIII. CONCLUSION 

In the present paper we extended Mott's and 
Shklovskii's approach for ID variable range hopping con- 
ductivity to short range interacting electrons in charge 
density waves and Luttinger liquids using an instanton 
approach. Following the recent paper^— , we calculated 
the quantum creep of charges at zero temperature and 
the linear conductivity at finite temperatures for these 
systems. The main results of the paper are equations 
(I4.16P for zero temperature and (j4.23p for nonzero tem- 
peratures, respectively. 

In our approach the quantum effects are weak but es- 
sential for the transport. An applied electric field renders 
all classical ground states metastable, and the motion 
between ground states along the electric field occurs via 
quantum tunneling. We were interested in the effect of 
the collective pinning by weak impurities, which always 
localizes a ID system with repulsive short range interac- 
tions on a sufficiently large length scale. The presence of 
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APPENDIX A 

Here we recall several relations for the bosonic repre- 
sentation of Id spinless fermions. In a one-dimensional 
system of short range interacting fermions, the excita- 
tions can be understood as density fluctuations of bosonic 
nature. Following Ref. 0, we define the compressibility 
as the derivative of the particle density with respect to 
the chemical potential: k = dp/dpi. Correspondingly, 
the elastic energy density is given by (Sp)"^ / (2k). 
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Using expression (|2.ip for the long-wavelength devia- 
tion of the density 6p — dx^/ir from po, we can rewrite 
the elastic energy in the form of the second term in the 
Hamiltonian (|2.7p if 



■nhv V 

, I.e. — 

K K nhK 



— . (A.l) 



The last combination is usually denoted by v^. 

As both pressure and chemical potential depend only 
on the ratio N/V, one can derive a relation between the 
conventional isothermal compressibility^^ 



]_dV 
VdP 



= P 



-2^ 



K 

72 



(A.2) 



and the above defined k. We should note that he com- 
pressibility in Ref. |20|was defined as p^{d^/dp) which is 
the inverse of the conventional thermodynamic isother- 
mal compressibility k. 

The kinetic energy is relevant for the description of dy- 
namical (quantum) phenomena. In a less rigorous way we 
may argue that the first term in Hamiltonian (j2.7p has its 
origin in a sum of single-particle kinetic energies p^/(2m), 
where p is the momentum conjugate to the displacement, 
which is proportional to ip, and m is the effective mass. 
Then, using commutator (|2.3p and the mentioned pro- 
portionality relations, we obtain a relation between the 
coefficient in the Hamiltonian and fermionic quantities 



Kv 



TThpo 



(A.3) 



This quantity is often referred to as vj. 

Thus, the parameters of the bosonized Hamiltonian, 
the stiffness K and the phason velocity, can be expressed 
in terms of the fermion system parameters'^: 



PqK 

m 




(A.4) 



In the (7-ology approach, the compressibility and the ef- 
fective mass are found as 



2n 



2lTpp 



TTh{2TTVF +54+52) 



2nvF +94 — 92 

(A.5) 

with the Fermi momentum pf — T^hpQ. Repulsion (K < 
1) usually corresponds to a smaller compressibility and a 
larger sound (phason) velocity. 



APPENDIX B 

In this appendix, we justify our choice of an Ohmic 
dissipation in Eq. (|3.19p and argue that it captures the 
generic features of more general, and possibly material 
specific, dissipative mechanisms. Coupling of the Id elec- 
tron systems to one- or higher dimensional phonon de- 
grees of freedom is the obvious candidate for the realiza- 
tion of dissipation. This mechanism is extensively dis- 
cussed in the literature on Mott variable range hopping. 



seeai for a review. The action for a bath with arbitrary 
dispersion relation and coupling to the electron system 
is given by 



•S'dis 



/3/2 



dT 



Ztt 



2"^"(^) [^Q2+^a ) Ca(T) + 



(B.l) 



Here, the Cq are bosonic bath degrees of freedom with 
eigenfunctions $q((7), eigenfrequencies Wq., a coupling 
U{q) to the electronic degrees of freedom, and 1/ilph 
is the volume element for one eigenstate in the a- 
representation. After integrating out the bath degrees 
of freedom, one finds the dissipative term 



Sdis 

IT 



-T 



^^{q)K{q,q';u;Mq') (B.2) 



with the kernel 



K{q,q;ujn) 



U{q)U{q') ^ <f>„(-g)$„(-g') 



2VL 



ph 



E 



. (B.3) 



The precise structure of the kernel K depends on the 
dispersion relation oja, the wave functions <I>„(g), and 
the coupling U{q) of the dissipative bath. For the in- 
stanton calculation in section IV we have used a generic 
Ohmic dissipation K{q,q';u!n) — r]6{q + q')\LL!n\ giving 
rise to a term rjLxhiLy in the instanton action. For a 
material specific choice of parameters one may not find 
this Ohmic dissipation, but instead a general dependence 
riF{Lx, Ly). In the following, we will argue that (i) for 
a weak coupling to the bath, one finds only a logarith- 
mic dependence of the instanton action on the external 
electric field and hence no change of the exponential field 
dependence in Eq |4.16l and that (ii) S'dis generally has a 
negative second derivative with respect to Ly. 

For a weak coupling to the bath, the dominant process 
will be the hop of a kink accompanied by the emission of 
one excitation of the bath. This process is described by 
a first order Taylor expansion 

The zero order term does not make any contribution to 
the imaginary part of the partition function Zi. In the 
spirit of such a Taylor expansion, one typically finds that 
the preexponential factor for Mott variable range hopping 
conductivity is proportional to the square of the electron 
bath coupling and in our notation proportional to 77. Ac- 
cording to this expansion, the instanton action contains 
a term In F{Lx, Ly) and depends via L^, Ly only in a 
logarithmic fashion on the external electric field. 

To illustrate our claim (ii) that generally ^ ''^g^T'^"'' < 
0, we consider the specific example of a local coupling of 
the Id electron system to Id acoustic phonons, i.e. the 
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general index a corresponds to momentum g, U{q) oc and hence a negative second derivative with respect to 
I/qQ^ (here Uq is the dimensionless electron-phonon cou- time. 
pUng constant), uja = "fph?- For an instanton with 
Lx ^ ~^^v oiic finds 



— cxconst.- -2-^ (B.5) 
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